Renormalization of the dephasing by zero point fluctuations 
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We study the role of zero-point-fluctuations (ZPF) in dephasing at low temperature. Unlike the 
Caldeira-Leggett model where the interaction is with an homogeneous fluctuating field of force, here 
we consider the effect of short range scattering by localized bath modes. We find that in presence 
of ZPF the inelastic cross-section gets renormalized. Thus indirectly ZPF might contribute to the 
dephasing at low temperature. 



I. INTRODUCTION 

The coherent motion of a quantum mechanical particle 
in a fluctuating environment is endangered by decoher- 
ence due to inelastic scattering events. The tempera- 
ture dependence of the resulting "dephasing" effect has 
been worked out in numerous studiesii^. During the last 
decade a controversy has emerged in the mesoscopic liter- 
ature regarding the role of zero-point-fluctuations (ZPF) 
in the theory of low temperature dephasing. The con- 
troversy was sparked by the experiment of Ref4> where 
a saturation of the dephasing rate in the limit of zero 
temperature has been reported, and consequently ZPF 
induced dephasing has been suggested as an explanation 5 
and debated^i^. 

Possibly one can insist that ZPF lead to "T = 0" de- 
phasing for a Brownian particle that interacts with an 
Ohmic Caldeira-Leggett (CL) bathiS where the fluctu- 
ations of the environment consist of long wavelength 
(q = 0) modes. Howeve r 11 ^ 2 , in metallic environment 
the effective fluctuations are characterized by a finite 
correlation distance, and hence consist of modes with 
wavenumbers q that range up to the Fermi momentum. 
It was largely accepted^ that if the interactions are 
short range, such that the fluctuations are characterize 
by a finite correlation distance, then the effect of the 
ZPF would be to renormalize the scattering cross-section 
and the mass of the particle. It turns out that in the 
cases of physical interest, and in particular for the pro- 
totype model of Refs i 15 i 16 , this renormalization effect is 
non-diverging in the zero temperature limit: both mass 
renormalization^ and dephasing 1? for a single particle 
in the presence of a dirty metal environment have been 
studied. Essentially the same formulation as in Refs^^ 
arise also in the more complicated many body treatment 
of the dephasing problem^. 

Though it is not diverging in problems of physical in- 
terest, still the renormalization effect due to ZPF might 
be significant in the actual analysis. The simplest possi- 
bility is to have an overall suppression of both elastic and 
inelastic scattering via a Debye- Waller factoid. But we 
would like to explore the more exotic possibility of having 
a distinct enhancement "factor" for the inelastic effect. 
It is therefore desired to have at least one exactly solvable 
model for dephasing due to short range scattering with 
environmental modes, that can be contrasted with the 



opposite CL limit where all the mode have q — 0. The 
objective of the present paper is to present such a model: 
In the proposed model (see Fig. 1) the environment con- 
sists of infinitely many localized fluctuating modes with 
(say) Ohmic spectral function, while the interaction is 
of short range and described by d (x) as in "s-scattering" . 
This should be contrasted with the long range interaction 
of the CL model which is linear in x. 

The outline of this paper is as follows: In Section 2 
we define a model for a localized bath that induces both 
zero point fluctuations (ZPF) and thermal fluctuations 
(TRF). In Sections 3-4 explicit expression for its scatter- 
ing matrix are derived following^. In Section 5-6 low 
temperatures are considered, where the TRF are treated 
as a small perturbation. The example and the numeri- 
cal analysis in Sections 7-8 establish that for weak TRF 
the effect of the ZPF background can be taken into ac- 
count by defining a renormalized intensity of the TRF. 
Accordingly ZPF may contribute to the dephasing at low 
temperatures, though not directly. 



II. THE MODEL 

The Hamiltonian of the particle plus the local bath is 

2m 



2 

n = ^- + s(x)J2c a Q a + J2 



n a 0J o 



(1) 



The index n a = 0, 1, 2, 3... may indicate the state of the a 
oscillator, or optionally n a — 0, 1 may indicate the state 
of a two level ("spin") entity, as in our numerics. From 
now on we use the notation 



Q = E< 



(2) 



Assuming an incident particle with kinetic energy we 
divide the oscillators into two groups: those with tu a < ek 
and those with cu a > e& respectively. We further assume 
low temperatures such that all the oscillators in the lat- 
ter group are in the ground state. Note that the particle 
has enough energy to induce real (non-virtual) excita- 
tion of any of the TRF oscillators. Hence we can write 
schematically: 



Q 



c s 



E - 

aeZPF 



E 

qgTRF 



(3) 
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where c s represents a static scatterer. The particle is 
affected by the fluctuations of Q. Assuming that the 
bath is prepared in the state n = m the fluctuations are 
characterized by the non-symmetrized power spectrum 



S(u>)= J2 \Qn m \ 2 27r5(u~(E n -E m )) (4) 

In the next section we show how the ZPF oscillators can 
be eliminated, such that the interaction is characterized 
by a dressed interaction matrix Q. Accordingly we define 
the effective power spectrum as 



S*(l>)= ]T \Q nm \ 2 2n5(LJ-(E n -E m )) (5) 
and the effective "size" of the elastic scatterer as 

C e ff = Qm,m (6) 

In the following sections we explain how to define Q and 
how to make the exact calculation of the elastic scat- 
tering amplitude T, and of the inelastic scattering cross 
section p inolastic (see Fig. 1). Then we discuss whether the 
results can be deduced from the effective values of c eff 
and S e{! (w). 

III. THE SCATTERING STATES 

Outside of the scattering region the total energy of the 
system (particle plus bath) is 



e fe + E n = e fc + 



n a uj a 



(7) 



We look for scattering states that satisfy the equation 

=£|tf) 

Open (propagating) channels are those for which > 
after scattering. Otherwise the channels are closed 
(evanescent). The channels are labeled as 



n = (no,n) = (no,n ZPP ,n T Rp) 
= (no, 711,712,713, ...,n a , ...) 



(8) 



where n = L,R for left/right, and 

^ZPF 1 ^TRF 3>I*6 collec- 
tive indexes for the two group of scatterers. We define 



k n = ^/2m(£ — E n ) for n € open (9) 
a n = \] — 2m(£ — E n ) for n G closed (10) 



later we use the notations 



v n = k n /m 
u n = a n /m 



(11) 
(12) 



and define diagonal matrices v = diag{w„} and u = 
diag{w„}. The channel radial functions are written as 

R(r) = A n e~ iknr + B n e +iknr n e open (13) 
R(r) = C„e- Q " r n e closed (14) 

where r — \x\. The wavefunction can be written as 

*(r,n ,Q) = Y,Rn 0<n (r)x n (Q) (15) 



The matching equations are 

*(0, right, Q)-*(0, left, Q) = 



(16) 



— [*'(0, right, Q) + *'((), left, Q)} = Q*(0,Q) (17) 
zm 

The operator Q is represented by the matrix Q nrn that 
has the block structure 



Qnrn 



Qw Q 
Quv Q 



(18) 



The matching conditions lead to the following set of ma- 
trix equations 

A R + B R = A L + B L 

-iv(A R -B R + A^-BJ = 2Q vv (A l , + B L ) + 2Q vu C l , 
-u(C Ti +C^ = 2Q uv (A lj +B lj ) + 2Q uu C lj 

From here we get the matching equations that relate the 
ingoing and the outgoing amplitudes: 

A R + B R = A L +B L (19) 
A r -B r + A Tj -B Tj = i2(v)- 1 Q{A^+Bj (20) 

where the dressed interaction matrix is defined as 

1 



Q — Qw Qv 



(U + Quu) 



Q, 



(21) 



In the next section we deduce the S matrix from the 
above set of equations, and obtain explicit expressions 
for the elastic scattering amplitude and for the inelastic 
cross section. 



IV. THE S MATRIX 

The unitary description of the scattering in terms of 
ingoing and out going probability currents requires to de- 
fine the normalized ingoing and the outgoing amplitudes 
as A n — ^/v~n~A n and B n = ^/vn~B n ■ Consequently we de- 
fined a re-scaled version of the Q nm matrix as follows: 



M, 



lvl uv lvl uu 



(22) 
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We also define a corresponding reduced matrix 

1 



M vu - 



-M u 



(23) 



Using these notations the set of matching conditions can 
be expressed using a transfer matrix as follows: 



B R 

ia 



(24) 



The transfer 2N x 2N matrix can be written in block 
form as follows: 

T. \ fl-iM -iM 
iM 1 + iM 



The S matrix is defined via 



(25) 



At 



B 

and can be written in block form as 

Sr St 
St Sr 

A straightforward elimination gives 
_ T -i T _ + 

r ++ -r_ + rrir + _ t+_ti 

Now we can write expressions for Sr and for St using 
the M matrix. 

St 



t -i 



(26) 



(27) 



(28) 



1 + iM 
St-1 



l-iM-M 2 + iM 3 + 



(29) 
(30) 



T = [S T ], 



The elastic forward scattering amplitude is 

1 

1 + iM_ 

The total elastic scattering probability is 

Elastic = m 2 + |T-l| 2 

= 1 - 2 [R(T) - |T| 2 ] 

We observe that the inelastic scattering is isotropic and 
its line shape (per direction) is 

p(w)= \[ST]nm\ 2 2tt6(lj - (E n -E m )) (33) 

with the measure cLj/(2tt). The total inelastic cross sec- 
tion is obtained by integration 



(31) 



(32) 



^inelastic 



- 2 J 



n(^m) 



1 



1 + 



(34) 



One can verify the Pi ne l as tic an< ^ ^elastic sum U P ^° um *y' 
which is essentially the "optical theorem" . 



V. PERTURBATION THEORY 

Since the temperature is low we treat the small effect 
of the TRF in leading order. We write 



Q = Q ZPF <g>l TE 
= Q° + SQ 



+ l z 



(35) 



where Q° is the sum over the ZPF coordinates including 
the static scatterer c s , while 5Q is the sum over the TRF 
coordinates. For the reduced Q matrix we get: 



Q = Q° vv + SQv 



-0° 



(36) 



u 



' uu 



U 



O 



u 



■> U U 



0° 



We assume that all the "important" open modes are well 
above the the evanescent threshold. This means that a 
single TRF transition is not enough to push the scattered 
particle into an evanescent mode. Accordingly 5Q UV and 
SQ VU are not included. With the same spirit we fur- 
ther assume that the TRF transitions hardly affect the 
evanescent velocity, hence 



^ZPF j^TRF 



(37) 



When calculating the matrix element Qnm the 
second term constitutes a sum over sequences 
Q^ v ...(5Q) v i.^...Q^ m . In order to have a non 
zero term, the TRF-oscillators of the v state should 
remain in the same state as in the n state, while one 
ZPF oscillator of the v state has to be excited. Similar 
observation applies to the states of the oscillators of the 
\x state. The TRF transitions are induced by SQ during 
the evanescent motion of the particle. Accordingly we 
deduce that 



S = 



+ 1 




ZPF 
UU 



ZPF 
UV 



Q 1 



u + Ql 

which can be written schematically as follows 

Q = c l TRF + A Q TRF (38) 

where the effective clastic scattering amplitude, and the 
scaling factor of the inelastic effect are 



co 



A 



ZPF 
VU 



ZPF 

VU 



u + Q 
1 

u + Q 



ZPF 
UU 



Q 



ZPF 

UV 



ZPF 
UU 



Q 



ZPF 
UV 



(39) 
(40) 



With an appropriate counter term we can make cq = 0. 
More interestingly we see that the effective TRF are char- 
acterized by the dressed power spectrum 



(Ao) 2 %) 



(41) 
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VI. THE DRESSED BORN APPROXIMATION 

The first order ( "Born" ) approximation relates the in- 
elastic line shape to the power spectrum of the fluctua- 
tions. We use the term "dressed Born approximation" 
in order to indicate that we use first order perturbation 
theory with respect to the TRF, while the ZPF including 
the static scatterer are treated to infinite order. Within 
this framework the leading order expression for the S 
matrix, using Eqs. (|38H40)) . is 

1 



1 + i{c /v £ ) + i\ Q M TRF 
= T l~tT 2 \ a M TRF + ... 

where the elastic forward scattering amplitude is 
% 



1 + i(co/v e ) 

Consequently we get for the inelastic scattering 
1 



|T | 4 (A ) 2 S{w) 



(42) 



(43) 



(44) 



Since we had assumed that the change in the kinetic en- 
ergy of the particle due to TRF inelastic scattering is 
relatively small, one can take v e - u ~ v e . 



VII. THE SIMPLEST EXAMPLE 

Consider a particle with velocity v e , that collides with 
a 'bath' that consists of an elastic scatterer c s , and a 
single two level TRF scatterer c T whose excitation energy 
is uj T (<ti efc). The interaction matrix is 



Q = Q 



C S C T 



(45) 



which we substitute in M « (l/w e ) 2 Q, so as to get 
St = (1 + iAi)^ 1 ■ In order to avoid crowded expressions 
we set the units such that v e = 1, and write 



1 



(1 + ic s ) 2 +c? 



-iOr 



— ICt 

1 + ic s 



(46) 



[Note again that in order to restore the units each c in 
the above expression should be replaced by c/v e ]. From 
here it follows that 



Pino 



(l-c 2 + i/ TRF ) 2 +4c 2 



[scaled] (47) 



= c 2 characterizes the intensity of the TRF, 
c s . One observes that for strong TRF the 



where u TBjF 
and Co = 

inelastic effect is suppressed and we get mainly elastic 
back reflection. But in the regime of interest, of weak 
TRF, the inelastic scattering is proportional to i/trf and 
agree with Eq.flU]) where \T \ 2 = 1/(1+Cq) and A = 1. 



Next we complicate the 'bath' by adding a single ZPF 
scatterer c z whose excitation energy is uj z (> e^). The 
possible values of the mode index are n = (0, 0) = 1, and 
n = (0, 1) = 2, and n = (1, 0) = 3, and n = (1, 1) = 4. 
The ZPF scatterer is assumed to be in the ground state 
(to = 1), and hence only the first two modes are open. 
The interaction matrix is 



Q 



I Cg C T C z 

c s c z 

Cz Cs c T 

\ c z c T Cs/ 



If we did not have the TRF oscillator, it would be a 2 x 2 
matrix: 



izrr | c s C 

. Cz Cs 



(48) 



If we ignored the ZPF oscillator, we would get Eq. (|4"B")) . 
But using Eq. ([2~Tj) we get the dressed interaction matrix: 



(U 3 +C S ) (W4+C S ) 



/ (u 3 +c s ) -c T 

I -C T (U4 + Csi 



with U3 = y/\ek ~ u 2 
sequently from c off = 



Cs 



| and u 4 = yj\tk- 
Qi,i we get: 

(u 3 +c s ) U ZPF 



and from u etf i 
A = 



(u 3 +c s ) (m 4 + 
Q2,i | 2 we get z/ eH 

^ZPF 



where 



{u 3 +c s ) (w 4 +c s ) - V T 



Con- 



(49) 



(50) 



Optionally we can get for 
Q the approximated result of Eq. (|3"5)) , which treats the 
TRF coupling in leading order. This treatment assumes 
that in the vicinity of the energy shell V\ w i>2 = v e , and 
u 3 m U4 = u z . The parameters Cq and Aq are calculated 



using Eqs. pHOl) with Q 
c z , leading to 



ZPF 
uu 



CO 



A 



c s 



(u z + Cs) 

^ZPF 

(u z + c s ) 2 



and - Q17 = 

(51) 
(52) 



In this simple example the dependence of Co and Ao on 
z/ ZPF is linear. But once we have more than one ZPF 
scatterer (as in the numerical example of the next sec- 
tion) the relation is no longer linear. It might be also of 
interest to solve the first equation c = for c s , and sub- 
stitute the result into the second equation. The outcome 
of this procedure is illustrated in Fig. 3. 

The calculation of the S matrix proceed in the same 
way as in the single TRF case, with the effective interac- 
tion matrix (no approximation involved): 



Q = 



(A-l)(u 3 +Cs) 
Ac-p Cs 



Ac T 

(A-l)(U4+Cs). 



(53) 
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Setting Di « i)2 e v e and U3 w U4 = u z as before, we label 
both diagonal terms as c off . Still we are not making any 
approximation with regard to the intensities v TRF and 
i^zpf, so as to get essentially exact results: 



k 2 + (c 8ff ) 2 R 2 



K 2 -(c ea ) 2 + A 2 i/ TRF ]2 + 4^( Ceff )2 



and the generalization of Eq. (j47| 
2 



-J |f | 2 |T| 2 A 2 ^TRF 



(54) 



(55) 



where |T| 2 is Eq.(|5"l"|l without the A 2 ^trf term. For weak 
TRF intensity, using A w Ao and T w 7o one obtains the 
dressed Born approximation Eq. ([44|) . One observes that 
the presence of the factor A has two implications: one is 
to enhance the inelastic scattering for weak TRF, while 
the other is to limit the range over which the weak TRF 
approximation applies. 



high temperatures (large !/ TRF ) this renormalization effect 
fades away and we get A = 1. See e.g. Eg. ([5(7)1 . The 
crossover involves a wild variation of A (see Fig. 4), which 
implies that < p in(! i a3 ti C < 1/2 goes through the whole 
range of possible values (Fig. 5). 

It is important to point out that if the fluctuations 
had continuous (rather than discrete) power spectrum, 
the above described intermediate wild variation would 
be smoothed away. Thus in realistic circumstances we 
expect that also in the presence of ZPF the qualitative 
dependence of p iaclaBtic on ^ TRF would be smooth, though 
renormalized by Ao at the limit of low temperatures. 

For the numerical study we consider a bath that con- 
sists of two level scatterers. The energy splitting of the a 
scatterer is u a and the interaction is described by the 
operator Q a =(° J). The strength of the interaction 
with the bath is characterized by the intensity of the 
fluctuations as obtained by integrating over their power- 
spectrum S(ui). Consequently we distinguish between 
the intensity of the ZPF and the intensity of the TRF: 



VIII. DISCUSSION, EXPECTATIONS, AND 
NUMERICAL DEMONSTRATION 



The analysis in the present paper is focused primar- 
ily on the low temperature scattering, due to weak TRF, 
where the dressed Born approximation of Section 6 ap- 
plies. Still in order to get the "big picture" we consider 
below the full range of ^ TRF values. We first highlight 
some qualitative observations that are based on the anal- 
ysis of the simple examples of the previous section, and 
then proceed with a numerical demonstration that in- 
volves a larger bath of scatterers. 

From the Born approximation we deduce that for 
weak TRF the inelastic cross section p illolaBtic is pro- 
portional to i/ THF . For strong TRF it drops down as 
implied e.g. by the simplest example Eq. ([4Tj) . The 
maximum p lnela . stlc = 1/2 is attained for the intermedi- 
ate value f trf = c o+l- A- priory we could not expect 
a larger inelastic effect because the elastic cross-section 
|T| 2 + \T— 1| 2 is bounded from below by the minimum 
value 50%. We can interpret the condition for attain- 
ing minimum elastic cross section using a Fabry-Perrot 
double barrier picture: The elastic scattering and the in- 
elastic scattering are like two barriers separated by an 
infinitesimal distance. The strongest interference effect 
is expected when the two barriers are comparable. 

The suppression of the inelastic effect for strong TRF 
is a generic effect: it becomes almost obvious if we con- 
sider the scattering of a particle from a fluctuating region 
in a three dimensional space. In the latter context strong 
fluctuations would repel the particle from the scattering 
region, hence making inelastic excitations within the ex- 
cluded volume less likely. So the strongest inelastic effect 
is experienced for intermediate values of v TRF . 

The inclusion of ZPF into the model renormalizes v TRF . 
The enhancement factor A is larger than unity (but finite) 
in the Born approximation limit, but if we go to very 



^zpf — 7 . C„ 
aeZPF 



E 

ezPi 

E ■ 

a£TRF 



(56) 
(57) 



The effective intensity of the thermal fluctuations is simi- 
larly defined and accordingly calculated from the dressed 
interaction matrix: 



E 



Or. 



(58) 



Given a set of N ZPF-scatterers with couplings c z , 
and a static scatterer c s , we calculate Co (which deter- 
mines Tq) and Ao as a function of v ZFF = -/V|c z | 2 . See 
Fig. 6. Then, for various values of v Z pf we calculate the 
exact results for v a[[ and plastic versus u TRF . See Fig. 7. 
One expects that for weak TRF the effective intensity v Bti 
would be proportional to v TRF , namely 



Kft ~ (Ao) 2 ^TRF 

Furthermore our perturbative scheme implies that 



^inelastic 



-2 m 4 

vi 



(59) 



(60) 



In order to test the quality of the latter approximation we 



re-plot the results for p u 



versus \T\ 



See Fig. 8. 



The numerical results confirm our qualitative expecta- 
tions, and are in agreements with the analysis of the 
simple example of the previous section. In particular one 
observes that the presence of ZPF has two implications: 
one is to enhance the inelastic scattering for weak TRF, 
while the other is to limit the range over which the weak 
TRF approximation applies. 



6 



IX. SUMMARY 

In the Caldeira-Leggett model the effect of the envi- 
ronment is characterized by a friction coefficient r\ and 
by a temperature T. But more generall y 9 ' 11 ' 1 ^ 18 it has 
been emphasized that the proper way to characterize the 
environment is by its form factor S(q,u>). The form fac- 
tor contains information on both the temporal and the 
spatial aspects of the fluctuations, and in particular one 
can extract from it not only T and r/, but also the spa- 
tial correlations. The general formula for the rate of 
dephasing 9,18 involves a dqdio integral over S(q,u>), and 
for short range interactions simply reflect the rate of in- 
elastic events. 

We find for our model system that the inelastic scat- 
tering cross-section p inclaatic is enhanced in the presence 
of ZPF, and accordingly ZPF may contribute to the de- 
phasing at low temperatures, though indirectly. This 
might come as a surprise since in HefJ^ it has been ar- 
gued that both elastic and inelastic scattering are sup- 
pressed by ZPF by the same Debye- Waller factor (DWF). 
A closer look reveals the difference between the two mod- 
els involved. In RefJ^ one considers the scattering of a 
particle (x) from a vibrating scatterer (Q), where the 
interaction is S(x—Q). Accordingly the particle experi- 
ences a fluctuating field U(x) = §{Q — x) and S(q,u) is 
the Fourier transform of (e~ l<?< 2W tfiQi )} ; which is sup- 
pressed by the DWF e - ^ q K In our model the interac- 
tion is 5{x)Q. Accordingly S(q,u>) is the Fourier trans- 
form of (Q(i)Q(0)), and, within the framework of the 
conventional Born approximation, there is no DWF in- 
volved: in our model adding high frequency components 
to the fluctuating field has no implication on the low fre- 
quency behavior of S(q,u>). 

The renormalization factor of the inelastic effect (A) 
in our dressed Born approximation comes from higher 



orders of perturbation theory with respect to the ZPF, 
while the TRF are treated in leading order. The renor- 
malization factor A multiplies the power spectrum S(u>) 
that describes the thermal fluctuations. The power spec- 
trum itself does not involve a DWF. The A renormal- 
ization of the inelastic scattering comes "on top" of the 
expected renormalization of the potential floor and of the 
inertial mass, which are familiar from the solution of the 
Polaron problem. In our scattering theory framework the 
expected renormalization of the potential floor can be de- 
duced from the ZPF induced offset in the effective 'size' 
of the elastic scatterer (c s ), while the renormalization of 
the mass comes from the associated energy dependence 
of the forward scattering amplitude (T). 

One can construct an extended bath that consists of 
an homogeneously distributed set of "s-scatterers" , as de- 
scribed inii. This would allow the modeling of a fluctu- 
ating environment of physical interest (say a Dirty metal 
environment) with the desired S(q,uj). In such physical 
circumstances we expect renormalization of (i) the po- 
tential floor; (ii) the inertial mass; and (iii) the effective 
thermal fluctuations. Our results imply that these renor- 
malization effects are non-divergent if the fluctuations 
are characterized by short range spatial correlations, but 
still they might modify the low temperature dependence 
of the dephasing effect. 
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Fig.l: Schematic diagram of the model system. The scattered 
wave of a particle that collides from the right with a thermal "s- 
scatterer" consists of forward elastic scattering with amplitude T, 
backward elastic scattering with amplitude T— 1, and isotropic in- 
elastic scattering with probability Pi no i aBt ic- Our purpose is to find 
the dependence of T and Pi nc i ast ic on the intensity of the low tem- 
perature thermal fluctuations (TRF), with arbitrarily large back- 
ground of zero point fluctuations (ZPF). 



/ 

8 




e 


















> 

time 





Fig. 2: Diagrams that describe the time evolution of kinetic energy 
can be used in order to illustrate terms the scattering calculation. 
The dotted line represents a contribution to the elastic cross section 
due to (virtual) scattering by ZPF modes. The solid lines represent 
contributions to the first order inelastic cross-section, where the 
intensity of the TRF is regarded as the small parameter. 
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Fig. 3: The renormalization factor for inelastic effect A ZPF is plot- 
ted as a function of vzpf for the simple model of Section 7, using 
Eg. 1521 with fixed static scatterer c s = (solid red curve) and with 
adaptive static scatterer such that co = (dashed blue curve) . The 
other parameters for this and for the next figures are u z = 0.96 and 
u>t = 0.03 and e k = 0.6. 
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Fig. 4: Plots of c e ff (upper panel) and A (lower panel) versus v T rf 
for the simple model of Section 7 with the same parameters as in 
Fig. 3. The red solid curves are for u ZPF = 0.5 and the blue dashed 
curves are for v Z pf = 5. 
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Fig. 5: The transmission |T| 2 (upper panel) and the inelastic cross 
section Pi no i as tic (lower panel) versus i/ TRF for the simple model of 
Section 7 with the same parameters as in Figs. 3-4. 
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Fig. 6: Plots of co and Ao versus Vzpf- For sake of comparison 
we also plot c G ft and A = (vaff/vTRF) 1 / 2 f° r * wo non-zero values of 
Ctrf. Here and in the next figure we consider a bath that consists 
of 7 TRF scatterers with u) a ~ 0.0003, and 4 ZPF scatterers with 
u> a ~ 0.96. There is no static scatterer (cs = 0). The kinetic energy 
of the incident particle is = 0.6. 
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Fig. 7: Plots of i/ B fi and Pmciautic versus Utrf for the same bath as 
in the previous figure. 
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Fig. 8: The inelastic cross section Pi no iastic versus the scaled inten- 
sity of the thermal fluctuations |T| 4 i/ c ff, using the data points of 
Fig. 7. 



